Non-perturbative interpretation of the Bloch vector's path 
beyond rotating wave approximation 
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The Bloch vector's path of a two- level system exposed to a monochromatic field exhibits, in the 
regime of strong coupling, complex corkscrew trajectories. By considering the infinitesimal evolution 
of the two- level system when the field is treated as a classical object, we show that the Bloch vector's 
rotation speed oscillates between zero and twice the rotation speed predicted by the rotating wave 
approximation. Cusps appear when the rotation speed vanishes. We prove analytically that in 
correspondence to cusps the curvature of the Bloch vector's path diverges. On the other hand, 
numerical data show that the curvature is very large even for a quantum field in the deep quantum 
regime with mean number of photons n < 1. We finally compute numerically the typical error size 
in a quantum gate when the terms beyond rotating wave approximation are neglected. 

PACS numbers: 42.50.Pq, 03.67.-a 
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I. INTRODUCTION 

The dynamics of two-level systems in an external field 
has been thoroughly studied for years. This problem ap- 
pears in several physical systems and is significant for the 
realization of quantum gates [2[ . The model consists of 
the Hamiltonian of a two-level system, coupled to a clas- 
sical field or to a quantum harmonic oscillator. There are 
no general methods for solving analytically these models 
without any approximations, so several approaches have 
been used to study their behavior. The best known is 
the Rotating Wave Approximation (RWA), that consists 
in neglecting the effects of the rapidly rotating terms [3| . 

For an atom residing in a resonant cavity the frequency 
fl of the Rabi oscillations between the two relevant states 
of the atom is typically 10 -6 of the atomic frequency uo a 
and of the cavity frequency cj, so that the RWA yields a 
good description of the system [4]. On the other hand, 
in circuit quantum electrodynamics (cQED) where 
superconducting qubits play the role of artificial atoms, 
one can enter the so-called ultrastrong coupling regime in 
which the ratio fl/u > 0.1 0-Q (at resonance, uj a = oj). 
In this regime, effects beyond RWA should be taken into 
account, both for the dynamics of pure states jlCH15j and 
for dissipative dynamics |16l-[l8j. 

The consideration of high Rabi frequencies is natural 
in quantum information theory, since high speed opera- 
tions are needed to perform a large number of quantum 
gates within the decoherence time scale, an unavoidable 
request for fault tolerant quantum computation p], 
Therefore, a deeper understanding of the effects beyond 
RWA is relevant for the prospects of quantum computa- 
tion. 

In this paper, we explain, in a non-perturbative man- 
ner, the non-trivial Bloch vector's path of a two-level 
system exposed to a monochromatic field. By consid- 
ering the infinitesimal evolution of the state vector, we 



determine the temporal dependence of the rotation axis 
and frequency. In particular, we show analytically that, 
when the field is treated as a classical object, the ro- 
tation frequency oscillates between zero and twice the 
value predicted by the RWA. The vanishing of the rota- 
tion frequency is associated with cusps in the Bloch vec- 
tor's path. Correspondingly, the curvature of the Bloch 
vector's trajectory diverges. On the other hand, we show 
numerically that the curvature can take very large val- 
ues even for a quantum field in the deep quantum regime 
with mean number of photons n < 1 . We finally show nu- 
merically that the size of the errors induced in a quantum 
gate by the terms beyond RWA scales as 

The paper is organized as follows. The model for a two- 
level system in a classical field is introduced and studied 
numerically in Sec.|TlJ while the analytical interpretation 
of the Bloch vector's paths is reported in Sec. HIB Nu- 
merical results for the quantum field model are shown in 
Sec. HV] where the errors introduced by the terms beyond 
RWA are also investigated. We finish with concluding re- 
marks in Sec. [Vj 



II. THE CLASSICAL FIELD MODEL 

We consider the following time-dependent Hamiltonian 
H(t), describing the interaction of a two- level system 
with a classical monochromatic field (we set h — 1): 

H(t)=H +Hi(t), 

J ff = e |0)(0|+ ei |l)(l|, (1) 

Hjit) = 2ficos(wt)(|0)<l| + |1)(0|), 

where u> is the frequency of the field and f2 the (Rabi) 
frequency of the field-induced oscillations between the 
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two levels |0) and |1) [19j. The time evolution of the 
two-level state vector \i/j(t)} = C (t)\0) + Ci(t)\l) is, in 
the interaction picture, governed by the equations 



iC (t) = n 

zCi(t) = n 



i(u) — UJa)t _|_ g — i(u)+UJa)t 



e i(uJ+U] a )t _|_ e -i{uJ-UJ a )t 



Ci(t), 
C (t), 



(2) 



where Co> a = ei — eo (ei > eo). 

The terms depending on uo + uo a oscillate very rapidly 
and are neglected in the rotating wave approximation. 
In this paper, we will explore the effects of these terms 
beyond RWA on the evolution of the two-level state vec- 
tor \ip(t)). For the sake of simplicity, we set the detuning 
A = uo — uo a = 0, as the treatment for A ^ would be 
essentially identical to that for A = 0. Finally, it is con- 
venient to normalize time in units of the Rabi frequency, 
that is, we set O = 1. Hence, we obtain 



iC = (1 + e- 2 ^ 1 ) Ci, 



ICi = (1 + e 



2iout\ 



(3) 



C . 



The RWA approximation is valid when uo ^> Q (in our 
units, when uo ^> 1). 

A convenient geometric picture of the evolution of the 
state vector is provided by the Bloch ball representa- 
tion p], Q , with the Bloch coordinates defined as 



X — 2Co r Cl r + 2CoiCli, 

Y = 2Co r Cii — 2CoiCi r , 
, Z = 2Co r + 2Cq { — 1, 



(4) 



where Ckr and Cki denote the real and imaginary parts 
of Ck (k = 0,1). The normalization constraint |Co| 2 + 
|Ci| 2 = 1 implies that the motion of the Bloch vector 
R = (X, Y, Z) takes place on the unit (Bloch) sphere 

x 2 + r 2 + z 2 = i. 

The most interesting feature the Bloch vector's path 
is the presence of cusps, shown in Fig. [TJ The corkscrew 
trajectories shown in this figure are obtained from numer- 
ical integration of the differential equations (|3]), starting 
from the initial condition Co(0) = 1, Ci(0) = (north 
pole of the Bloch sphere). The plots of Fig. [1] show, from 
top to bottom, the Bloch vector's trajectory and its pro- 
jections on the yZ, XZ, and XY planes; from left to 
right, uj = 2.5, 5, 20, and RWA approximation (uo — >> oo). 
While within RWA the trajectory is a circle, correspond- 
ing to Rabi oscillations between states |0) and |1) (north 
and south pole of the Bloch sphere, respectively), cusps 
appear beyond RWA. We will interpret in a simple but 
exact manner these numerical results in Sec. IIII1 



III. INTERPRETATION OF THE BLOCH 
VECTOR'S PATH 

We consider the infinitesimal evolution of the state vec- 
tor: 



\ijj(t + dt)) = U(t,t + dt)\i/)(t)), 



(5) 



where the unitary operator U(t,t + dt) represents an in- 
finitesimal rotation of the Bloch vector through an an- 
gle dO about the axis directed along the unit vector 
n = (nj,ny,nz): 



dO 

U(t,t + dt) = ( J-iy n-cr 



(6) 



where / is the identity operator and a — (ax, cry, (?z), &k 
(k = X, y, Z) being the Pauli operators. In the {|0), |1)} 
basis, Eq. ([5j) reads as follows: 



C (t + dt) 
d(t + dt) 

1-inzf -(n Y + in x )f 

(riY-inx)f 1 + inzf 



Co(t), 
Ci(t) 



(7) 



Since 



\ij>(t+dt)) = \ij>(t)) + \i>(t))dt, 

with \ip(i)) obtained from Eq. (j3j), we have 
'C (t + di) 



(8) 



d(t + dt) 
1 



-i (1 + e 
1 



-2iujt\ 



c (t), 

Ci(t) 



(9) 



From comparison between Eqs. © and (j9j) we obtain 
dO 



< 



n x - 



[1 + cos(2a;t)]dt, 



dO 



nyy = sin(2o;t)dt, 



(10) 



K n z = 0- 

From the first two equations of this system we have 
d6 



(11) 



Therefore, the rotation speed vanishes when cos(cjt) 
0, that is, for 



t = tk 



(2fe + l)7T 

2^ ' 



(k = 0, 1, ...). 



(12) 
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FIG. 1. Bloch vector's path (top plots) and its projections on, from the second to the fourth row, the YZ, XZ, and XY planes. 
From left to right: uj = 2.5, 5, 20, and RWA (cj —> oo). Time evolution is followed, starting from the north pole, up to time tv 
(within RWA, at such time the Bloch vector returns, for the first time, back to the north pole). 



After insertion of Eq. ([TT]) into Eq. ([TO]) , we obtain 
' nx = | cos(c<;£)|, 

< riy = $m(ojt) sgn [cos(o;t)], (13) 
, n z = 0, 

where sgn(x) = x/\x\ is the sign function. Therefore, at 
times tk the rotation axis changes instantaneously from 
(nj,ny, nz) = (0, 1, 0) to (0, —1, 0), that is, from the Y 
axis to the opposite direction — Y. Such discontinuity in 
the rotation axis is associated with a cusp in the Bloch 
vector's path. As discussed in appendix [Aj the curvature 



of the Bloch vector's path diverges at the cusps. 

The dependence of the rotation axis and speed on time 
can be visualized by means of the schematic drawing in 
Fig. [21 Within RWA, the rotation speed is given by the 
length #rwa = 20 = 2 (O = 1 in our units) of the seg- 
ment OA and the rotation axis is the X axis. On the 
other hand, when effects beyond RWA are taken into ac- 
count, the rotation speed at time t is given by the length 
= 4| cos(co>£)| of the segment OB, and the rotation axis 
is directed along OB. Therefore, the rotation speed os- 
cillates between = (at times tk) and = 4 = 2#rwa 
(at times tk = kir/uj, k = 0, 1, ...). 

In Fig. [3j we show the XY-plane projection of the 



FIG. 2. Schematic illustration of the rotation axis and speed 
within and beyond RWA. 



Bloch vector's evolution in the neighborhood of a cusp 
and compare it with the RWA approximation. It is in- 
teresting to remark that the distance between the two 
evolutions, shown at different times by dashed lines, is 
not simply given by the shortest distance between the 
exact trajectory and its RWA. This is due to the fact 
that, as pointed out above, while the RWA predicts a 
constant rotation velocity #rwa = 2, the exact evolution 
exhibits a variable rotation speed = 4| cos(co>£)|, which 
is faster than in the RWA far from the cusp and slows 
down up to — at the cusps. 

The above discussed non-trivial behavior of the Bloch 
vector's path implies, as shown in Fig. HJ the existence of 
plateaus in the dependence of the path length s(t) as a 
function of time. The plateaus are obtained around the 
times tk where the curvature of the Bloch vector's path 
diverges. The slope of s(t) oscillates between zero (at the 
cusps) and 2#rwa = 4. Note that the overall path length 
is larger than within RWA, since the Bloch vector moves 
along a path longer than within RWA (see Fig. [3j). The 
mean slope of s(t) approaches the RWA slope in the limit 

UJ — » 00. 

Finally, we point out that, while the curvature exhibits 
a singularity at the cusps, the functions X(i),Y(i), and 
Z(t) are regular at the times t = tk corresponding to such 
cusps. This is due to the fact that, as discussed above, 
the rotation speed tends to zero when t —> tk HO]. 



IV. THE QUANTUM FIELD MODEL 

The interaction between a two-level system and a sin- 
gle mode of the quantized electromagnetic field is de- 
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FIG. 3. Comparison between the XY projection of the Bloch 
vector's exact evolution (full curve) and its RWA (straight 
line) in the vicinity of a cusp, for uj = 5. The dashed lines 
link points corresponding to the exact and approximate Bloch 
vectors, R(t) and Rrwa(£), obtained after the same evolution 
time t (the different dashed lines correspond to different values 
oft). 
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FIG. 4. Plateaus in the path length s(t) as a function of time 
(full curve), for uj — 5. The straight dashed line corresponds 
to the RWA result, s(t) = 2t. 

scribed by the Hamiltonian (h = 1) Q 
H = H o + Hj , 

H ® = \ u aVz + u [a ] a + ^ , (14) 

Hi = Aa + (a 1 " + a), +AV_ (a 1 " + a), 

where a± = | (ax ± ^y) are the rising and lowering op- 
erators for the two- level system: <r + |0) = |1), cr + |l) = 0, 
<r_ |0) = 0, cr_ 1 1) = |0), the operators and a cre- 
ate and annihilate a photon: a^\n) p = \fn + l\n + l) p , 
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a\n) p = y/n\n — \n) p being the Fock state with n 
photons. For simplicity's sake, we consider the resonant 
case (lj = uo a ) and the coupling constant A G M. The 
RWA is obtained when we neglect the term cr+a\ which 
simultaneously excites the two-level system and creates 
a photon, and a_a, which de-excites the two- level sys- 
tem and annihilates a photon. In this limit, Hamiltonian 
([T4]) reduces to the Jaynes-Cummings Hamiltonian (see, 
for instance, Ref. [3]). Within RWA (Jaynes-Cummings 
limit), there are coherent Rabi oscillations between the 
two- level system-cavity states |0)|n + l) p and \l)\n — 1} P . 
The frequency of such oscillations is Q n = Xy/n. 

We study numerically the temporal evolution of the 
two-level system's state vector, when the initial state is 
|0o)|a)p? where |0o) = cos (|) |0)+sin (|) |1) is the initial 
state of the two- level system and \a) p = Xl^Lo c ^l n )p' 
with c n = exp(— - |a| 2 )-^=, is the coherent state of the 

field, with mean number of photons given by n = \a\ 2 and 
root mean square deviation in the photon number An = 
\J~Tt. The classical field limit is obtained for ft — >• oo, 
with constant Q, n = X^/E = Aa, that is, with A ex 1/a. 
We measure frequencies in units of fi^, namely we set 
Qn = 1, i.e., A = 1/a. 

We focus on the time evolution of the Bloch ball coor- 
dinates. In Fig. [5]we show the Bloch vector's trajectory 
and its projection on the YZ, XZ, and XY planes, both 
in the deep quantum regime a = 1 (first two columns 
in Fig. [5]) and closer to the classical field limit (a = 5, 
corresponding to n = 25 photons, in the third column 
of Fig. [5]). For small values of a the Bloch vector's path 
crucially depends on the initial state |0o) of the two- 
level system, as we can see from the two cases shown 
in Fig. [5j 9 = (north pole of the Bloch sphere) and 
6 = 7r (south pole) (such strong dependence does not 
appear for large values of a). In all instances, id = 5, 
so that the effects beyond RWA are pronounced. The 
plots for a = 5 are quite close to the classical field plots 
(fourth column in Fig.[5j). On the other hand, cusps, even 
though smoothed, are very pronounced even in the deep 
quantum regime a = 1 . In Table [I] we show the value of 
the curvature k at times t\ = ^ and ti = where n 
diverges for the classical field model (see appendix |A|) . 
Note that the obtained values are much larger than the 
value k = 1 corresponding to RWA in the classical field 
model. 





first peak 


second peak 


a = 1,0 = 
a = 1, = 7r 
a = 5,0 = 


2.6 x 10 4 
12 

1.5 x 10 5 


1.8 x 10 2 

8.5 
2.8 x 10 3 



TABLE I. Value of the curvature k for the first two peaks, at 
times t\ — and ti — f^j where k diverges when the field 
is treated as a classical object. The three lines of this table 
correspond to the first three columns of Fig. [5] so that uo = 5. 

We finally estimate the error introduced by the terms 



beyond RWA in a typical quantum gate. Since the time 
elapsed during an elementary quantum gate is of the or- 
der of the inverse of the Rabi frequency, we compute the 
Euclidean square distance on the Bloch ball between the 
exact and the RWA evolution, 

||R - Rrwa|| = 

(15) 

y/(X - Xrwa) 2 + (Y- Yrwa) 2 + (Z- Zrwa) 2 

as a function of time, up to t = r = tt/Q. The root mean 
square 



5 = \j\[ *H R W " R *waW|| 2 (16) 

indicates the typical size of the errors introduced in a 
quantum gate by the terms beyond RWA, once such 
terms are neglected. Our numerical data shown in Fig. [6] 
show that 8 ex 1/uo for any value of a. Moreover, the 
values obtained for a quantum treatment of the field ap- 
proach quite quickly the expectation of the classical field 
model, and already for a = 5 we obtain values of S close 
to the classical values. 



V. CONCLUSIONS 

In this paper, we have explained features beyond RWA 
of the Bloch vector's path of a two-level system exposed 
to a monochromatic field, including the oscillation of the 
rotation frequency and axis and the presence of cusps in 
the trajectory. Our exact and non-perturbative analysis 
is based on the infinitesimal evolution of the two-level 
system state vector. 

The above non-trivial features of the Bloch vector's 
path should be taken into account, together with the ef- 
fects of noise, in any implementation of quantum gates, 
where error rates smaller than 10 -4 are requested for 
fault-tolerant quantum computation. In particular, since 
cusps are clearly seen even for the quantum field in the 
deep quantum regime with a small mean number of pho- 
tons n < 1, they should be relevant and observable 
in circuit QED experiments. The numerical results of 
Fig. [6] show that neglecting the terms beyond RWA one 
introduces an error of the order of the inverse of the ra- 
tio between the field and the Rabi frequency. There- 
fore, to obtain error rates smaller than 10 -4 , effects be- 
yond RWA should be necessarily taken into account when 
Q/uj > 10 -4 , a limit largely overtaken in the ultrastrong 
in which the ratio Q/uj > 0.1. Finally, the existence of 
cusps due to the presence of two frequencies, in this case 
u — u a and id + id a , appears as a rather general phe- 
nomenon, so that we expect cusps to appear in protocols 
such as the Raman transition [21]. This remains to be 
analyzed in future [221 ] . 
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FIG. 5. Bloch vector's path (top plots) and its projections on, from the second to the fourth row, the YZ, XZ, and XY planes. 
From left to right: a = 1,0 = 0; a = 1, # = 7r; a = 5, # = 0; and classical field, = 0. In all instances, uj = 5 and time evolution 
is followed up to time tt. 



Appendix A: Curvature of the Bloch vector's path 



Given a point R = (X, Y, Z) of the Bloch vector's tra- 
jectory, we compute its velocity as 



V = R = Rx6>n, 



(Al) 



with n = (nx^TiY-nz) direction of the rotation axis and 
rotation frequency. After substitution of Eqs. ([TT]) and 



([T3]) into (jAip , we obtain the components of vector V: 
Vx = — 4Zsin(co>£) cos(cjt), 

V Y =4Zcos 2 (cjt), (A2) 
^ Vz = 4cos(cj£) [Xsin(cjt) - Y cos(ujt)}. 
We then compute the tangent vector 

dn v 



t = 



ds 



(A3) 
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is given by 



FIG. 6. Root mean square distance between the exact evolu- 
tion and the RWA, for, from bottom to top, a = 1 (circles), 
2 (squares), 5 (diamonds) and classical field (triangles). The 
straight full line shows the slope. 



where s is the arc length along the Bloch vector's path 
and 



s = |V| = 4| cos(ut)\y/Z 2 + [Xsm(ut) - Y cos(cjt)] 2 . 

(A4) 

Finally, the curvature of the Bloch vector's trajectory 



|R- st| 



(A5) 



with the components (Ax,Ay,Az) of the acceleration 
A = R given by 



' A x = -4Zcjcos(2u;t), 

Ay = -4Zu)8Ul(2u)t), 

< A z =4£j[Xco8(2u)t) + Ysm(2u)t)], 



(A6) 



and 



5 = 4 ( — | cos(out)\j ^Z 2 + [Xsm(ut) — Y cos(ut)} 2 
4 1 cos(cj£)| 



^Z 2 + [Mpj - Ycos(ujt)} 2 



x - [(X 2 - Y 2 ) sin(2cjt) + XY cos(2a;t)]. 

2 (A7) 
As shown in Fig. [7J the curvature diverges at the times 
tk of Eq. ([T2]h corresponding to the cusps in the Bloch 
vector's path, i.e., when cos(oot) = 0. 
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